Nuclear parton distributions at next to leading order. 
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We perform a next to leading order QCD global analysis of nuclear deep inelastic scattering and 
Drell-Yan data using the convolution approach to parameterize nuclear parton densities. We find 
both a significant improvement in the agreement with data compared to previous extractions, and 
substantial differences in the scale dependence of nuclear effects compared to leading order analyses. 
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I. INTRODUCTION 

Ever since the discovery that quark and gluons in 
bound nucleons show momentum distributions notice- 
ably different to those measured in free or less bound 
nucleons two decades ago, the precise determination 
of nuclear parton densities has awakened growing atten- 
tion, driving both increasingly precise and comprehen- 
sive nuclear structure functions measurements 2], and 
more refined theoretical understanding of the underlying 
physics. 

The precise knowledge of nuclear parton distribution 
functions (nPDF) is not only required for a deeper under- 
standing of the mechanisms associated to nuclear binding 
from a QCD improved parton model perspective, but it 
is also the starting point for the analyses of a wide variety 
of future and ongoing high energy physics experiments, 
such as heavy ion collisions at RHIC Q, proton- nucleus 
collisions to be performed at LHC |J| , or neutrino- nucleus 
interactions in long baseline neutrino experiments , for 
example. Consequently, the accuracy of nPDF is rapidly 
evolving into a key issue in many areas of particle physics. 

From the point of view of perturbative QCD, the ex- 
traction of nPDF can be done in close analogy to what 
is done for free nucleons: they are considered as non per- 
turbative inputs, to be inferred from data, whose rela- 
tion to the measured observables and energy scale de- 
pendence are computed order by order in perturbation 
theory. Although one can not discard larger higher twist 
power corrections than in the case of free nucleons, or 
even some sort of nuclear recombination effect, factoriza- 
tion and universality of nPDF are expected to hold in a 
very good approximation, over a wide kinematical range 
in the present experiments. 

At variance with nucleon PDF, which driven by the 
demand of more and more precise predictions, have at- 
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tained in the last years an impressive degree of accuracy 
and refinement, nPDF extractions are in a considerably 
earlier stage of the development. Not only the num- 
ber and diversity of nuclear data is much more reduced, 
but the analyses are restricted to leading order accuracy 
(LO), with rather crude parameterizations of nuclear ef- 
fects which lead to poor /d.o.f values in global QCD 
fits to data 0,01. 

There are also some caveats inherent to the particu- 
lar approaches implemented so far, that define nPDF in 
terms of nucleon PDF and a multiplicative nuclear cor- 
rection factor at a given initial energy scale, from where 
they are evolved, fn the first place, nPDF defined in 
this way have their momentum fraction per nucleon re- 
stricted to be less or equal than unity, what excludes from 
the description, and also from the evolution, a portion of 
its natural kinematical range. Of course, unless one is 
specifically interested in that region, this approximation 
is expected to imply a minor correction, although one 
can not discard the excluded region as a source of new 
and interesting information. 

Much more problematic in these approaches is that 
the actual shape of the nuclear correction factor required 
to reproduce accurately the data implies very capri- 
cious functions, with a large number of parameters, and 
which in practice precludes the numerical computation of 
the scale dependence at next to leading order accuracy 
(NLO). The actual computation of structure functions 
and evolution equations at this order implies several con- 
volution integrals very difficult to evaluate unless Mellin 
transform techniques arc applied. 

In this paper we show that a much more convenient 
alternative to deal with nuclear effects is to define nPDF 
using a convolution approach. In such framework the free 
nucleon parton densities are convoluted with very simple 
weight functions that parameterize nuclear effects. The 
convolution method naturally takes into account the ac- 
tual range of nucleon momentum fractions, allows via 
Mellin transform techniques a straightforward numerical 
evaluation of the NLO scale dependence, lead to extraor- 
dinarily accurate nPDF with relatively few parameters, 
and finally, allow to interpret the nuclear modifications 
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in terms of a very simple mechanism of rebalance of mo- 
mentum fractions between the distributions. The success 
of the convolution approach comes from the fact that 
the momentum fraction dependence of nuclear effects is 
strongly correlated to that of partons in free nucleons, as 
shown in re-scaling models 8] , feature which is explicitly 
included by the convolution. 

As an example of feasibility of the convolution ap- 
proach we obtain for the first time a full NLO ex- 
traction of nPDF from a large number of nuclear DIS 
El III El m and DreU Yan data jl^. We also asses the 
differences between the LO and NLO extractions, find- 
ing that although the quality of fits to present data is 
comparable in both approximations, there are important 
differences between the dependence of the nuclear cor- 
rection factors using either LO or NLO extractions. This, 
for example, questions the use of LO factors with NLO 
nucleon parton densities to generate nPDF in NLO com- 
putations, as it is common practice. 

In the following section we summarize the main moti- 
vation and features of the convolution approach. In the 
third section we present the details and outcome of the 
NLO and LO nPDF extractions, comparing our results 
with previous LO analyses. In the last section we dis- 
cuss the differences between LO and NLO extractions of 
nPDF, computing the one hadron production cross sec- 
tion with them as and example, and present our conclu- 
sions. 

II. NPDF 

The description of DIS processes ofi^ nuclear targets, 
eA e' X, is customarily done in terms of the hard scale 
Q^, defined as minus the virtuality of the exchanged pho- 
ton, and a scaling variable xa, analogue to the Bjorken 
variable used in DIS off nucleons, 

^ -q\ XA ^ ^r^. (1) 

2pA-q 

respectively. Here pA is the target nucleus momen- 
tum and consequently, xa is kinematically restricted to 
< < 1, just as the standard Bjorken variable. 
Alternatively, one can define another scaling variable 
xn = Ax A, where A is the mass number of the nucleus. 
Under the assumption that the nucleus momentum pA 
is evenly distributed between that of the constituent nu- 
cleons pn — Pa/A, this variable resembles the Bjorken 
variable corresponding to the scattering off free nucleons, 
xjv = /{2pffq). However in nuclear scattering context 
of course it spans the interval < x^ < A, by definition, 
and reflecting the fact that a parton may in principle 
carry more than the average nucleon momentum. 

When discussing nPDF the usual approach is to pro- 
pose a very simple relation between the parton distribu- 
tion of a proton bound in the nucleus, f.f, and those for 
free protons fi, 

ifixN^Ql) = R^{xN,Ql A., Z) MXN,QI), (2) 



in terms of a multiplicative nuclear correction factor 
Ri{xj^, Q^, A, Z), specific for a given nucleus {A, Z), par- 
ton flavor i, and initial energy scale Qq. This description 
is convenient since the ratio Ri(x]\j,Q^,A, Z) compares 
directly the parton densities with and without nuclear 
effects, and is closely related to the most usual nuclear 
DIS observables, which are the ratios between the nu- 
clear and deuterium structure functions. Indeed, we will 
use this picture to simplify the presentation of the out- 
put of our analysis. However, this is not the best suited 
way to parameterize the effects in the intermediate steps 
of the analysis, nor the best alternative for higher order 
numerical computations, as it was mentioned previously. 

Since fi{xN,Q'^) are deflned for < a; at < 1, nPDF 
deflned as in Eq.Q are a priori restricted to that same 
range at the initial scale. In most analyses the evolu- 
tion equations are also constrained to < a; at < 1, loos- 
ing the possibility of analyzing nuclear effects beyond 
xn > 1 at any other scale. Notice that even if nPDF 
were constrained to be zero for xat > 1 at a given scale, 
the complete evolution equations would produce non-zero 
values for other scales in the range of xn neglected by 
the approach, leading to a possible violation of momen- 
tum conservation. From another side, the direct param- 
eterization of the ratios Ri{xN , A, Z) require many 
parameters and even XTv-dependencies whose Mellin mo- 
ments can not be put into a closed expression. 

A much more adequate alternative to Eq. (0) is to relate 
nPDF to standard PDF by means of a convolution 

ff{xN, Ql) = r ^W,{y, A, Z) h(^,Ql). (3) 

A™ V y 

where the weight function Wi{y,A, Z) parameterize now 
the nuclear effects and can be thought as the effective 
density of nucleons within the nucleus carrying a frac- 
tion y/A of its longitudinal momentum. Besides allow- 
ing the full kinematical range for nPDF, this kind of 
approach has been shown to lead to remarkably good 
parameterizations of nuclear effects with very few pa- 
rameters, and with a smooth A-dependence For ex- 
ample, neglecting nuclear effects, the effective nucleon 
density is just Wi{y,A,Z) — A6{1 — y), and simple 
shifts in the momentum fraction carried by nucleons 
Wi{y, A, Z) — A5[l — y — e) have been shown to repro- 
duce in a very good approximation many features of the 
known nuclear modifications to structure functions. 

The convolution approach also allows to perform the 
evolution in Mellin space, which is much more convenient 
numerically, and almost mandatory for NLO accuracy. 
Defining yA = ^ Eq.® reads 

f,^{AxA,Ql)^ f ^wMyA,AZ)M^,Ql) 
JxA yA yA 

(4) 

and in Mellin space, i.e. taking moments in both sides of 
Eq.Q, 

f^{N, Ql) = W^{N, A, Z) h{N, Ql) (5) 
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where 



and similarly for /i(iV, Qq) with A = 1, and 



(6) 



W^{N,A,Z)= dxx^-^W^{Ax,A,Z). (7) 
Jo 

Notice that the moments are defined in terms of the cor- 
rect scaling variable xa-, and that the evolution of the 
nPDF moments ffiN, Qq) can be managed with a stan- 
dard evolution code in Mellin space. 

Nuclear structure functions are defined as the average 
of the proper combination of the bound proton and neu- 
tron structure functions as 

AF^ix.Q'') = ZF^^^ix,Q^) + {A-Z)F^'/^ix,Q^), 

(8) 

where both bound nucleon structure functions can be 
written in terms of the corresponding nuclear parton dis- 
tributions in Mellin space as 



(9) 



2tt 



The first term in the right hand side of Eq.® corre- 
sponds to the leading-order contribution while the second 
represents the next-to-leading order correction. Expres- 
sions for the NLO quark and gluon coefficients Cq^\N) 
and dg^\N) can be found in In our analysis we 

consider only 3 active flavors and neglect the contribu- 
tion from heavy quarks to the structure function. 

Considering that DIS data does not allow to determine 
all combinations of flavors, and that nuclear effects may 
be expected to be isospin invariant in a first approxima- 
tion, it seems to be reasonable to introduce only three in- 
dependent Wi{y, A, Z): one for the valence distributions, 
another for the light sea, and the last for the gluons. In 
this way 



fo) 



11) 
11 



) 

d^{N,Ql) 



= W^{N,A,Z) u,{N,Ql) 

= W,{N,A,Z) d,{N,Ql) 

= W,{N,A,Z) u{N,Ql) 

= W,{N,A,Z) diN,Ql) 

= WgiN,A,Z) giN,Ql) 



(10) 



Notice that in many proton parton densities, and in par- 
ticular in GRV98 JLl], the one will be used in the follow- 
ing, the strange and heavy quark densities vanish at the 
low initial scale and there is no need to introduce addi- 
tional weight functions for the latter. The parton distri- 
butions for the bound neutron can be obtained from the 
ones above by isospin symmetry. 



An interesting feature in this approach is that since 
non-singlet combinations of parton densities, as those 
for valence quarks, evolve independently, and both 
ff{N, Q^) and /i(-/V, Q"^) obey the same evolution equa- 
tions, the moments of the weight function Wi{N, A, Z) 
are scale independent . This is not the case for sea or 
gluon weight functions, since they are affected by singlet 
evolution resulting in a scale dependence due to quark 
and gluon mixing and therefore needed to be defined at 
a particular scale. 

Charge, baryon number and momentum conservation 
imply constraints to be satisfied by the weight functions: 

W,iN = l,A,Z)^l (11) 
u^i2,Ql)+d^i2,Ql) 
+2u^{2, Ql) + 2d^{2, Ql) + 5^(2, Ql) = 1 (12) 

The best results are obtained using for the valence dis- 
tributions weight functions like 



W,{y,A,Z) = 
A [a^ (5(1 



(13) 



y) + (1 - a„) (5(1 - - y)] 



A) 



A) 



where the first two terms may be interpreted as reduc- 
tion in the parent nucleon longitudinal momentum frac- 
tion, and in spite of their simplicity reproduce accurately 
EMC and Fermi motion effects. Indeed, with just these 
three parameters the fit reproduces fairly well the large 
xn data, however it fails to account for antishadowing 
effects at intermediate xjq, where valence distributions 
still dominate. 

In order to include antishadowing effects we add the 
third term in Eq. (|13|l . which induces a small enhance- 
ment of the distributions with > and a mild xn 
dependence given by the parameters and Notice 
that the convolution integral dilutes the xn dependence 
producing an effect similar to those predicted by parton 
recombination models I* would be pointless to try 
to extract from the data the nuclear modification to va- 
lence distributions in the low xn region, since this region 
is clearly dominated by sea and gluon densities, however 
the first three terms in Ea. H13|) violate charge conserva- 
tion. In order to remedy this situation we include the 
forth term, similar to the third but where charge conser- 
vation forces < 0. In this case the weight function 
can not be interpreted as a probability density but as 
the result of the mechanism that compensate re-scaling 
and recombination effects. The parameters as and [3s 
are taken to be the parameters used for the sea quark 
densities, fixed by the low xjsi behavior of the data, while 
both and Us are fixed by momentum and charge con- 
servation, Eas. lfTT)) and ifT^ . 

For sea quarks and gluons, the argument is just the op- 
posite: the fit to data is not sensitive to any nuclear mod- 
ification at high or intermediate a; at, so the best choice 
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for sea weight functions is found to be 
Ws{y.A, Z) = A5{1 -y) + ^^ (f (l - f )''^ (14) 

where the first term gives the gluon distribution in free 
protons and the eulerian function (specifically the pa- 
rameter as) drives the low Xn shadowing (a^ is found to 
be negative). The three parameters, a^, tts and fig are 
strongly constrained by the structure functions ratios at 
small xm and the Drell-Yan ratios. Ng is just the normal- 
ization of the eulerian function, Ns = B{as + 2,(3s + !)■ 
Similarly, for gluons 

W.,(,,A,Z)=A.(l-,) + ^(l)"^(l-l)'^'.(15) 

Here, Og is strongly constrained by the Q^-dependent 
F^"^ jF^ data, whereas the exponents are taken to be 
the same as for the sea distributions and Ng — Ng. No 
significant improvement is found assigning independent 
parameters for them. Since Uy and are constrained 
by momentum conservation in Eas. Hll(l and (|12(l . there 
are 9 independent parameters e^, e„', a^, a^, Pv^s, cts, Ps, 
and ttg for each nucleus. Since no data is available on 
different isotopes of the same nucleus, in the following 
we drop the dependence on Z . The A dependence of all 
the parameters can be written 

= 7, + X,A^^ (16) 

The full parameterization of both the A and xn depen- 
dence of nuclear effects imply in principle quite a lot of 
parameters (27), however the mild A dependence found 
in some of them allow to reduce their number. For exam- 
ple, Us , ay and Uy can be taken as constant in A loosing 
just a few units in but eliminating 6 parameters. 

III. RESULTS 

In the following section we present results from the 
LO and NLO fits to nuclear data using the convolution 
approach. The data analysed include the most recent 
NMC and SLAC-E139 DIS structure functions ratios to 
deuterium and carbon, selecting those measurements cor- 
responding to Q > 1 GeV, and also E772 Drell-Yan dilep- 
ton cross sections from proton nucleus collisions, as listed 
in Table 1, rendering a total number of 420 data points. 

The parameters were fixed minimizing the function 
defined in the customarily way as 

where af^ stands the measured observable, cr*'' the cor- 
responding NLO (LO) estimate, the statistical and 
systematic errors added in quadrature, and the sum runs 
over all the data points i included in the fit. No artificial 
weight has been given to any particular subset of data. 



TABLE I: Nuclear data included in the fit. 



Measurement 


Collaboration 


Refs. 


7^ data 




NMC 


[13] 


18 




SLAC-E139 


[14j 


18 


Fp/Fi' 


SLAC-E139 


[14j 


17 


F?/Fi> 


NMC 


[13] 


18 




SLAC-E139 


[14j 


7 


Ff/Fi' 


SLAC-E139 


[14j 


17 


Fp/Ff 


NMC 


m 


18 




SLAC-E139 


[14] 


7 




SLAC-E139 


[14] 


23 


Fi'/F^ 


SLAC-E139 


[14] 


7 


F^^/F^ 


SLAC-E139 


[14J 


18 


F.p/F.f 


NMC 


[15] 


15 


Ff/F^ 


NMC 


[15] 


15 


Fp/F^ 


NMC 


[15] 


15 


Fi'VFf 


NMC 


[15] 


15 


fP/f§ 


NMC 


[15] 


15 


Fp/F§ 


NMC 


[16] 


145 


"dyI"dy 


E772 


[17] 


9 


Ca in 
"DYl"DY 


E772 


[17] 


9 


„Fe iD 

"dy/"dy 


E772 


[17j 


9 


W iD 
"DYI"DY 


E772 


fl7] 


9 


Total 






420 



In order to compute the structure function and Drell- 
Yan cross-section for Deuterium, for which we neglect 
the nuclear corrections, and as a basis to construct the 
observables for the different nucleus, we use the parton 
distributions in the free proton as provided by the GRV98 
analysis ^Jl- Consequently, we fix the initial scale to 
Ql ^ 0.4 GeV2 (0.26 GeV^) and the QCD scale (for 

five flavors) A^^^^^ = 167 MeV (132 MeV) to NLO (LO) 
accuracy. The total obtained is 300.15 units for the 
NLO fit and 316.35 for LO one. 

The comparison between the data on the ratios of dif- 
ferent nuclear structure functions to deuterium and those 
computed with NLO nPDF is shown in Figure 1 (the LO 
prediction one is almost indistinguishable). The solid 
line corresponds to the result of the fit computed at the 
value of each experimental point, whereas the ex- 
trapolation to small Xn has been performed using the 
corresponding of the smallest- a; at point. For heavy 
nuclei, the low a; at is mainly constrained by DIS ratios 
to carbon, as shown in Figure 2. The structure function 
ratios are useful to determine mainly the valence quark 
distributions, while Drell Yan data, shown in Figure 3 be- 
comes crucial in order to fix sea quark distributions. The 
gluon distribution enters the structure function at NLO 
or through the scale dependence of the parton densities, 
making very difficult to obtain information about it in 
DIS experiments. The dependence oi F2"' / F2 , shown 
in Figure 5, is the most sensitive available observable to 
the gluon distribution. Nevertheless, it is worth pointing 
out that there is still a large uncertainty on this density 
and data from observables where the gluon distribution 



5 




FIG. 1: F^jF^ data. The lines interpolate the values ob- 
tained with the NLO nPDF set at the respective Q^, and ex- 
trapolate to low XN at the leftmost point. 




0.8 



Anmc 

Ca/C 



0.8 



0.8 






1 10"^ 10"' 1 



FIG. 2: The same as Fig. 1 but for F^/F? data 



enters at the lowest order, like in hadronic colliders, is 
needed to obtain a much better constraint. 

The regular A dependence of the parameters, as ob- 
served in Figure 5, helps to interpolate through regions 
where the data is scarce and also lead to reasonable ex- 
trapolations where there is not available. 




FIG. 3: Data on nuclear Drell "fan cross sections rates to 
deuterium and those computed with NLO nPDF. 



Noticeably, while some parameters show a clear depen- 
dence on the size of the nucleus, such as the shifts in the 
momentum fractions e and e' which drive nuclear effects 
at moderate and large x^, those related to the shape of 
the nucleus effective densities at small x^, such as a,,, 
as and = as are not strongly dependent on A. The 
well known A dependence of shadowing effects at small 
xiq is driven by the normalization of these effective den- 
sities as, and a^, and also by the large xn behavior of 
the densities fixed by the parameters (i^ and /3s, which 
control how much of the large Xj^ component of the PDF 
enters the convolution. 
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FIG. 4: Scale dependence of Fi"^ IF§ data and the outcome 
of NLO nPDF. 
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FIG. 5: A -dependence of the parameters. 
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FIG. 7: Ratios coming from different nPDF sets 



The resulting nPDF are shown in Figure 6, as ra- 
tios to free proton PDF, as defined in Eq.(|21l, for va- 
lence quarks, light sea and strange quarks, and gluons, 
at = 10 GeV^. The numerical computation of these 
ratios once the nPDF have been written and extracted 
within the convolution approach, is straightforward and 
allows a comparison with standard analyses, and other 
parton distribution functions. The ratios are also pro- 
vided (in a FORTRAN code) as grids in xm,Q'^ and A 
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FIG. 6: The ratios for the densities are computed at 10 GeV^ 



for practical purposes and can be obtained upon request 
from the authors. 

Similar results are found using nucleon parton densities 
other than GRV98. The similarity between modern par- 
ton densities guarantees that the nuclear ratios obtained 
with a given PDF set lead to reasonable nPDF when 
combined with another. Therefore, the parton distribu- 
tion in a proton of a nucleus A can be simply obtained 
by multiplying the nuclear ratios obtained in our anal- 
ysis by any modern set of proton PDF. Of course, this 
is true provided the distributions come from analyses at 
the same order in QCD, as we will discuss in the next 
section. 

It is worth mentioning that the agreement between nu- 
clear parton distributions and data is remarkably better 
in the case of convolution based parameterizations than 
the ones found with multiplicative parameterizations in 
previous (LO) analyses. Those analyses yield values 
around 630 in the case of Ref. (6| , and even larger values 
with the parameterization of Ref.0, for the same data 
set used in the present analysis. A detailed comparison 
with the parameterization of Ref.^] indicates that the 
of that set is basically due to a large contribution 
from the Q^-dependence of / F2 data. Therefore, 
one might expect significant differences between our LO 
gluon nuclear ratio and that of EKS98. 

In Figure 7 we show a comparison between the nuclear 
ratios for LO Uvaienco, u and g proton (in calcium) distri- 
butions from nDS, EKS98 and HKM sets, at = 2.25 
GeV^ and 100 GeV^. 

The main differences with previous nuclear parameter- 
izations are found in gluon densities, and to less extent in 
sea quarks, when comparing with ESK98, and in valence 
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and sea quarks in the case of HKM. Our gluon densities 
show comparatively less small xjq shadowing for heavy 
nuclei than ESK98, and very little antishadowing at in- 
termediate a:; AT, which is considerable in their distribu- 
tions. For nuclei lighter as C our gluons show only a tiny 
antishadowing effect. 

In order to study the sensitivity of different observables 
on the amount of shadowing in the nuclear gluon distri- 
bution, we have performed an alternative extraction of 
nPDF from the same data set but constraining the gluon 
density in heavy nuclei to show a stronger shadowing ef- 
fect at small xm- We provide the result in set called 
nDSg, constrained to satisfy i?^ " — 0.75 at xm = 0.001 
and = 5. The value of this analysis is around 
550, considerably larger than the unconstrained fit, and 
should be considered only as a mean to study variations 
on, mainly, the gluon nuclear distribution. An example 
for the gluon nuclear ratio in Au is shown in the last plot 
in Figure 6. 

Compared with HKM, our valence distributions show 
low xn shadowing, whereas in HKM there is none. In 
our parameterization the ratios for sea and gluon densi- 
ties approach unity as x^ grows, while with HKM distri- 
butions these ratios show a strong rise. The difference in 
the fits may be understood due to the fact that both the 
low xn region in valence densities, and the large xm be- 
havior of the sea distributions have little impact in DIS 
observables, and is only picked up by Drell-Yan yields, 
not included in HKM analysis. 



IV. NLO 

Although there are no significant differences between 
the total values obtained in LO and NLO fits, the cor- 
responding nPDF, and also the ratios to ordinary PDF, 
Ri{xj\r,Q^,A), indeed differ, specially at low Q^. In some 
cases, the differences are as large as the nuclear effects 
themselves. In Figure 8 the ratios between the NLO and 
LO extractions of R{xn, Q^, A) 



Kf{xN,A,Q'' 



R: 



NLO 



(xAr,Q2^A) 



Rl^'ixN.Q^A) 



(18) 



are shown as a function of a; at for various and different 
nuclei. 

The main differences between Rf (x n , Q'^ , A) and 
Ri^ {xn , A) are found in sea quark and gluon den- 
sities, most noticeably at low xjy and in heavy nuclei. 
These differences are correlated through the scale depen- 
dence with a similar behavior in the valence ratios at 
small xn- As one would expect, the differences are more 
significant at small and fades away as increases. 

The fact that LO analyses lead to fits comparable in 
accuracy to NLO analyses is due to the relatively mod- 
erated range in spanned by the data, and the absence 
in the data set of nuclear observables strongly dependent 
on the gluon distributions. This, of course neither im- 
ply that the differences in R{xntQ^,A) as obtained at 
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FIG. 8: Ratios between NLO and LO extractions of nuclear 
effects for different nuclei. 



LO and at NLO accuracy would be negligible, in fact be- 
yond LO the ratios are factorization scheme dependent, 
nor that there would not be significant NLO corrections 
in other observables. This is particularly true for those 
that rely on the gluon or sea quark densities and stress 
the importance of having NLO extractions of nPDF. 

As an example of an observable which is sensitive to 
NLO corrections, in Figure 9 we show the LO and NLO 
leading twist cross sections for the production of neutral 
pions in d-Au collisions as a function of the transverse 
momentum px of the final state particle. The plots cor- 
respond to nucleon center of mass energies ^/snn of 200 
GeV and a range in pseudo-rapidity of jryj < 0.18, com- 
puted using the code in adapted for nuclear beams. 
Pion fragmentation functions were taken from Ref . [l9| . 
The LO and NLO cross section are computed at two dif- 
ferent values for the factorization and renormalization 
scales = = Pt and fi^i — M_f — '^Pt- The differ- 
ences between both predictions give an estimate of the 
theoretical uncertainty in the fixed order calculation. As 
can be observed, there is a considerable reduction in the 
scale dependence of the cross section when the NLO cor- 
rections are taken into account. This feature is found 
in almost any infrared-safe observable in hadronic colli- 
sions, indicating that LO calculations can only provide a 
qualitative description. The inset in Figure 9 shows the 
Knlo factor, defined as 



K 



NLO 



..NLO 

'dAu 

LO ■ 
'^d.Au 



(19) 



i.e. the ratio between the NLO and LO cross section, 
computed at a given factorization and renormalization 
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FIG. 9: LO and NLO neutral pion production cross sections 
for d-Au collisions. 



scales, and using the corresponding nPDFs and fragmen- 
tation functions. Notice that the K^lo factor is not 
a physical quantity, actually it turns out to be strongly 
scale dependent, but provides an estimate of the size of 
the NLO corrections. In this case, it shows that the one- 
loop QCD corrections for pion production at RHIC are of 
the order of 50% or even larger, while the strong increase 
in the correction when the transverse momentum pT is 
smaller than ~ 2 GeV indicates we approach the limit of 
the region of validity of pQCD calculations. 

Finally, in Figure 10 we show the estimates for the 
d-Au cross sections to neutral pions but normalized to 
proton-proton cross sections, against the data reported 
by |20j, and not included in the present fit. With the ex- 
ception of the lowest pr data, which is in the borderline 
of the perturbative domain, the estimate agrees with the 
data well within the experimental uncertainties. Notice 
that even though the LO and NLO cross sections dif- 
fer substantially, the discrepancies almost cancel in the 
ratios, to proton, highlighting the perturbative stability 
and the consistency between the LO and the NLO extrac- 
tions of nPDF. To quantify the effect of nuclear shadow- 
ing in the gluon distribution on hadron production at 
RHIC we have computed the same observable using the 
alternative NLO set of nPDF (set nDSg) with a larger 
shadowing in the gluon distribution at small x^- The 
Figure shows a reduction in the d-Au cross section at 
small pt for the set with larger gluon shadowing, how- 
ever it seems unlikely to obtain a much smaller value 
for the ratio RdAU a.t pr ^ 2 GeV with reahstic nPDF. 
Notice that the cross sections receives contributions typi- 



FIG. 10: Neutral pion production nuclear rates. 



cally from xn > 0.01. Similar conclusions, with a slightly 
bigger reduction at small pt, are reached if the same ob- 
servable is computed for neutral pion production in the 
forward region (77 ^ 3). Any experimental finding on a 
stronger reduction of RdAU might certainly be consid- 
ered as an evidence of a new phenomena, at least beyond 
pQCD. 



V. CONCLUSIONS 

We have performed for the first time a full NLO QCD 
global analysis of nuclear DIS and Drell Yan data using 
a convolution approach to parameterize nPDF. We have 
found that this strategy not only leads to much more 
accurate nPDF but simplifies considerably the numeri- 
cal computation of QCD corrections at NLO. Although 
both LO and NLO nPDF reproduce the available data 
with comparable precision, they show non negligible dif- 
ferences which have to be taken into account when com- 
puting other observables. 
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TABLE II: Parameters of the NLO and LO nPDF 



Parameter 


1 


NLO 
A 


5 


7 


LO 
A 


5 




0.1984 


-0.0013 


0.0814 


0.2030 


-0.0014 


0.9510 




0.0346 


-0.0124 


0.9421 


0.0351 


-0.0133 


0.3657 




0.7546 


-0.6687 


-0.0473 


0.7251 


-0.6647 


-0.0583 




2.1412 


2.2690 


-0.0390 


2.1786 


2.5720 


-0.0439 




-0.0474 


0.3730 


1.1301 


19.925 


2.2760 


1.1463 


as 


-0.0135 


-0.0202 


0.2797 


-0.0179 


-0.0189 


0.2664 


Cts 


0.7980 


0.0814 


-0.8647 


1.0616 


0.0572 


-0.6277 


Ps 


-24.325 


7.3191 


1.1204 


-24.107 


7.3526 


0.4284 


ttg 


0.0565 


-0.0073 


0.4244 


0.0629 


-0.0076 


0.4285 


x'/d.o.f. 




299.91/393 






316.35/393 
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